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THE AVERAGE PHYSICAL PROPERTIES AND STAR FORMATION HISTORIES 
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ABSTRACT 

We investigate the average physical properties and star formation histories (SFHs) of the most 
UV- luminous star-forming galaxies at z ~ 3.7. Our results are based on the average spectral energy 
distributions (SEDs), constructed from stacked optical to infrared photometry, of a sample of the 
1,913 most UV-luminous star-forming galaxies found in 5.3 square degrees of the NOAO Deep Wide- 
Field Survey. We find that the shape of the average SED in the rest-optical and infrared is fairly 
constant with UV luminosity: i.e., more UV luminous galaxies are, on average, also more luminous 
at longer wavelengths. In the rest-UV, however, the spectral slope (3 (= dlogF\ / diogX; measured at 
0.13/im < A rcs t < 0.28/im) rises steeply with the median UV luminosity from —1.8 at L « L* to —1.2 
(L ks 4 — 5L*). We use population synthesis analyses to derive their average physical properties and 
find that: (1) L\jy, and thus star-formation rates (SFRs), scale closely with stellar mass such that more 
UV-luminous galaxies are more massive; (2) The median ages indicate that the stellar populations 
are relatively young (200-400 Myr) and show little correlation with UV luminosity; and (3) More UV- 
luminous galaxies are dustier than their less-luminous counterparts, such that Lk4 - 5L* galaxies 
are extincted up to A(1600) = 2 mag while L m L* galaxies have A(1600) = 0.7 — 1.5 mag. We argue 
that the average SFHs of UV-luminous galaxies are better described by models in which SFR increases 
with time in order to simultaneously reproduce the tight correlation between the UV-derived SFR 
and stellar mass, and their universally young ages. We demonstrate the potential of measurements of 
the SFR-Af» relation at multiple redshifts to discriminate between simple models of SFHs. Finally, 
we discuss the fate of these UV-brightest galaxies in the next 1 — 2 Gyr and their possible connection 
to the most massive galaxies at z ~ 2. 

Subject headings: cosmology: observations — dust extinction — galaxies: evolution — galaxies: 
formation — galaxies: high-redshift — galaxies: stellar content 



1. INTRODUCTION 

Constraining the main mode of star formation at high 
rcdshift is of fundamental importance to the theory of 
galaxy formation. Significant progress has been made 
in characterizing some of the general physical proper- 
ties of the star-forming galaxy population at high red- 
shift. Measurements have been made of the cluster- 
ing, spectral energy distribu tions (SEDs), and ki n emat- 
ics of these galaxies fe.g.. [Shaplev et al.l 120011. 120031: 
Papovich et al.l 120011: iPentericci et al.l 120071: ILee et al.l 
20061 120091 : IStark et all 120091 : IVanzella et al.l 120091 ). De- 
spite these efforts, major questions remain about the de- 
tails of how these galaxies form their stars. We do not 
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yet know whether their star formation history is domi- 
nated by a series of short but intense "bursts" or a long 
continuous accretion of gas converted into stars. Con- 
straints on the star formation history of galaxy popula- 
tions can bear on the availability of fuel for star forma- 
tion, and hence their gas accretion history. For example, 
bursts may point to merger- induced processes which de- 
liver galaxy-sized parcels of gas Ce.g.lMihos fc Hernquistl 
119961 iKolatt et al.lll999llSomerville et al.ll2001h . whereas 
continuous star formation might suggest a more sta- 
ble fu el delivery system , such as "cold-mod e accretion" 
(e.g., iKeres et al.l [20081: iDeke! et al.l [200l) Observa- 
tions at z = - 1 (e.g., iNoeske et aL]|2007t lElbaz et ah! 
I2007h and a limited set of observations at redshift out 
to z ~ 2 (|Genzel et al.l [2001 iDaddi et al.ll2007l ). how- 
ever, show some incons istencies wit h theoretical predic- 
tions at high redshift (Dave 2008). For example, the 
studies of faint star-forming galaxies at z > 3 seem to 
favor short duty cycles unlike theoret ical expectations 
(e.g., ILee et al.l 120091 : IStark et al.l[2009h . 

It is currently unclear whether such discord requires 
changes to the theory (e.g., implementing luminosity- 
dependent galaxy formation) or is the result of obser- 
vational biases inherent in our incomplete census of 
galaxies at high redshift. This is because most re- 
cent efforts have been focused on samples derived from 
the small-area, deep surveys (e.g., GOODS, HUDF: 
iSawicki fc Thompson! [2001 iBouwens et al.ll2007t ). which 
while excellent at identifying sub-L* galaxies, lack the 
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cosmic volume to contain significant samples of the 
most luminous and massive galaxies. While the UV- 
faint galaxies likel y dominate the cosmic star formation 
rate density (e.g., M adau et al.lll996t iReddv fc Steidell 
2009), their measured star-formation r ates, space den- 
sities and clu stering amplitudes (e.g.. iLee et al.l 120061 : 
iBouwens et al.ll2007l:lReddv fc Steidelll2009ft suggest that 
they are not likely to be the proge nitors of the most 
massive galaxies seen at later epochs dQuadri et al.ll2007t 
lElsner et all 120081 : iMarchesini et al.l 12009ft . In contrast , 
the UV-brightest galaxies at z > 3 may provide a more 
plausible progenitor population for the most massive 
galaxies (M sta r ^ a few x 10 11 M Q ), especially if they are 
able to sustain their observed star-formation rates for an 
extended duration (> 1 Gyr). 

In this paper, we define a large sample (1913) of the 
most UV luminous galaxies at z > 3, and investigate 
their average properties with the goal of determining 
their stellar masses, star formation histories, and evo- 
lutionary descendants. 

We use (fl,n A ,a 8 ,h wa ) = (0. 28,0.72,0.9,0 72). M ag- 
nitudes are given in AB system (O ke fe Gunnlll983l ) un- 
less noted otherwise. 

2. DATA AND GALAXY SAMPLE 

Our sample of 1,913 super-L* star-forming galax- 
ies at z rs 3.7 is selected from the 9.3 deg 2 Bootes 
field of th e NO AO Deep Wide -Field Survey (NDWFS, 
hereafter: Uannuzi fc Devi 11999ft . The NDWFS dataset 
used in this work inc ludes the optical By/RI imaging 
(jJannuzi fc Devi ll 9991 available from the NO AO Science 
Archive) taken with the MOSAIC camera, near-infrared 
JHK S taken with the NEWFIRM camera on the May- 
all 4m telescope (A. Gonzalez et al., in prep), and the 
Spitzer IRAC data from the Spitzer Deep Wid e-Field 
Survey (3.6, 4.5, 5.8, 8.0 ^m: lAshbv etall [20091 avail- 
able from the Spitzer Heritage Archive). Where avail- 
able, we also use [/-band ("[/-spec" filter) data taken 
with the Large Binocular Telescope (F. Bian et al., 
in prep) to check whether galaxies in our sample are 
detected. Because the optical data were taken over 
multiple pointings, the depth and seeing vary signifi- 
cantly over the entire Bootes field. In order to en- 
sure robust photometric measurements and candidate se- 
lection, we restrict our study to areas where the see- 
ing in the Bw, R and / bands has a Moffat pro- 
file with a full width at half maximum MFWHM< 
1.3", with a best-fit Moffat parameter (3 < 2.5. Us- 
ing these criteria, we exclude five subfields of the ND- 
WFS data (NDWFS J1426p3421, NDWFSJ1428p3456, 
NDWFS J1428p3531 , NDWFSJ1431p3346, and ND- 
WFS J1434p3456). The seeing and photometric depths 
of the optical data in all subfields are found on the sur- 
vey websit43- After excluding these subfields and the 
masked area, the effective area of our sample is «5.3 
degree 2 . 

Our sample is selected by applying a Lyman-break 
color-selection technique to the BwRI data. The 
Lyman-break technique is designed to identify a UV 
bright star-forming galaxy with a strong spectral break 
(at wavelength A < 1216 A) resulting from absorption by 
the intervening Lya forest. At redshifts 3.3 < z < 4.3, 

14 http : //www . noao . edu/noao/noaodeep/ 



this break lies between the Bw and R ba nds. Follow- 
ing n u merous studies in the l i terature (e.g ., Stcide fet al.l 
119991: iGiavalisco et al.l 12004 IBouwens et al.l 12007ft . we 
simulated the expected colors of typical star-forming 
galaxy SEDs at a range of redshifts and defined the fol- 
lowing color criteria: 

(B w - R) > 1.2 + 3.0 x (R - I) (B w -R)> 2.0 
(R - 1) > -0.3 S/N(R) > 3 S/N(I) > 7 (1) 

In order to only include objects with robust photomet- 
ric measurements, we excluded candidates lying within 
masked regions around bright stars (including diffrac- 
tion spikes around them), large galaxies, and the edge 
of each subfield image. We also excluded objects 
that are heavily blended with adjacent sources and 
stars based on the SEx tractor CLASS_STAR indices 
(|Bertin fc Arnoutsl fl996ft . Of 108 objects that satisfy 
the criterion CLASS.STAR > 0.95, 78 objects (72%) are 
brighter than I = 23.0 where the stellarity index can be 
measured reliably. Furthermore, we use optical-infrared 
colors (R—I vs. I— 3.6//m define a ste llar locus well sepa- 
rated from galaxies; see Figure 8b in iBrown et al.ll2007ft 
to further exclude stars whenever possible. We also vi- 
sually examined the LBT J7-band data and exclude any 
of candidates which are detected; the [/-band samples 
light below the Lyman continuum break (A = 912 A) 
at redshifts z > 3.4, and hence any candidate detected 
in U is unlikely to be at the redshift of interest. Very 
few (<1%) of the sources are identified in the U-band 
(with detections defined as 3er or higher, corresponding 
to U < 26 mag). After excluding these individual de- 
tections, we stacked the resulting sample and confirmed 
that the median candidate is also not detected in the 
(stacked) /7-band image (see Figure [4j. We note that 
less than 10% of our candidates are formally detected in 
the IRAC channel 1 and 2 (3.6/mi), and less than 1% 
in channel 3 and 4. While we do not exclude any ob- 
ject based on the IRAC photometry, this is reassuring 
because the relatively shallow depth of the Spitzer data 
of the Bootes field is unlikely to warrant detection of 
most high-redshift sources except for very massive ones 
(see later). The 5-sigma depths of the SDWFS data are 
22.74, 22.24, 20.39, 20.25 mag (AB) for IRA C channel 1, 
2, 3, and 4, respectively (jAshbv et al.ll2009ft . 

The final number of candidates is 1,913, correspond- 
ing to a surface density of ss 0.10 arcmin~ 2 over the 5.3 
degree 2 area. As shown in Figure Q] (top right), our sam- 
ple contains a significant number of very bright sources 
- for example, 236 sources brighter than Iab — 23.5; 
these provide ideal targets for detailed studies of galaxy 
properties such as their stellar populations, interstellar 
media, and kinematics. 

3. SPECTROSCOPIC SAMPLE 

In order to characterize the redshift distribution and 
contamination rates by low-redshift interlopers and 
Galactic stars, redshifts for a small sample of candidates 
were obtained using the Deep Imaging Multi-Object 
Spect rograph spectrograph (DEIMOS; iFaber et"aT1 
12003ft on the Keck II Telescope of the W. M. Keck Obser- 
vatory on U.T. 2010 May 15. Three masks, containing 
a total of 66 candidates, were observed with 2.0 — 2.5 
hours each using the standard DEIMOS setup with the 
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Fig. 1. — Left: The B\^RI color-color plane with the adopted color selection criteria for z > 3.3 candidates denoted by the black lines. 
The black dots represent a subset of all NDWFS galaxies. Colored symbols indicate the sources with spectroscopic redshifts, demonstrating 
the relatively clean selection of high-redshift targets. Dashed lines show the theoretical tracks of star-forming galaxies at z > 3.3 (z = 3.3 
at the bottom) for four reddening parameters (E(B — V)=0.0, 0.15, 0.30, 0.45). Right: The magnitude distribution of the full sample 
( "FULL" : solid) and spectroscopic sample ( "SPEC" : dashed) are shown (top) . The bottom panel shows the redshift distributions of the 
entire spectroscopic sample (dotted line) and the subset containing only the fainter galaxies with I > 22.5 (solid line). We adopt the latter 
as the representative redshift distribution for the full photometric sample. The expected redshift distribution of the full sample from our 
simulations is in good agreement with the observed one (dashed line). 



600 1/mm grating (Ablaze = 7500 A) centered on ~ 7600 
A. Mask positions were chosen to maximize the num- 
ber of primary targets (secure -Bjy-band dropout candi- 
dates; S/N > 8) and secondary targets (fainter candi- 
dates; 5 < S/N < 8) within the DEIMOS field of view. 
The spectra were reduced using the DEIMOS pipeline 
developed by the DEEP2 survey team at UC Berkeley 
(M. Cooper et al., in prep). While the details of all our 
spectroscopic observations will be presented elsewhere, 
we show typical ID spectra of our candidates in Fig- 
ure [H In addition to these observations, we compiled 
existing spectroscopic redshifts from the past Keck ob- 
servations by Spinrad, Soifer and collaborators as well 
as those obtained from the AGN and Galaxy Evolution 
Survey (AGES: C. Kochanek et al., in prep). 

In total, 78 sources from our sample have been ob- 
served spectroscopically. Of these, redshifts could not 
be determined for 12 sources due to optical faintness 
or a lack of useful lines. Two of the spectroscopically 
confirmed sources a re Syitzer MIPS 24/im sources at 
z=4.230 and 3.675 (|Desai et alJl20"08h . The former is 
a narrow-line AGN exhibiting CIV and Hell, and the 
latter is a broad-line QSO. Additionally, three sources 
have been identified as QSOs from our DEIMOS spec- 
tra (z=3.835, 4.230, 4.050), and one source is a Galactic 
star. The total number of confirmed interlopers (i.e., 
those other than star- forming galaxies at z ~ 3.7) is 5 
AGN and 1 star. Hence, the rate of contamination ranges 
from 8% (if all 12 unidentified sources are at high red- 
shift: 6/78) to 23% (if all 12 sources are at low redshift: 
(6+12)/78). In Figure [1] {bottom right), we show the two 
redshift distributions of our spectroscopic samples, one 
for the full sample including the MIPS sources, and the 
other showing only / > 22.5 galaxies without MIPS de- 
tections. We further note that the MIPS 24 /mi sources 
are fairly bright (/ < 22.5) compared to the majority of 



our candidates. Hence, we adopt the latter as the rep- 
resentative redshift distribution N(z), as the majority of 
our sample is fainter than / = 22.5. Thus, the redshift 
distribution is well described by the mean (z) = 3.7 and 
the standard deviation of 0.4. 

To further validate the adopted redshift distribution, 
we carried out photometric simulations to estimate the 
expected redshift distribution of our dropout sample 
based on the broad-band color selection. The simulation 
consists of 1) creating artificial star- forming galaxies over 
a range of redshifts, UV luminosity, spectral slope, mor- 
phology, and size, 2) inserting them in the real images, 
3) color-selecting in an identical manner in which the 
re al galaxies were select ed. We refer in t ereste d readers 
to iFerguson et alJ (|2004l ) and lLee et all (|2006| ) for more 
details. The expected redshift distribution is shown in 
Figure [T] (bottom right; dashed line), in good agreement 
with the observed one. 

4. MEASUREMENTS OF THE AVERAGE SPECTRAL 
ENERGY DISTRIBUTIONS 

We investigate the average SED properties of the UV- 
selected star-forming galaxies by stacking our imaging in 
the B W RIJHK S [S.6][A.5][5.8][8.0] passbands on the po- 
sitions of the 1,913 galaxies in our sample. For a galaxy 
at z — 3.7 - the median redshift of the sample - these 
bands probe the rest-frame wavelength range 0.1 — 2/tm. 
In order to investigate whether their UV brightness (a 
proxy for star formation rate) is correlated with other 
galaxy properties such as stellar mass or extinction, wc 
define six subsamples ordered by their J-band magni- 
tudes. In Table 1, we summarize the details of each sub- 
sample. Given the relatively narrow redshift distribution 
(z = 3.7 ±0.4), the /-band flux provides a fairly accurate 
representation of the observed UV luminosity measured 
at the rest-frame 1700 A. At a fixed /-band flux, the UV 
luminosity of a source at z — 3.7 will be higher (lower) 
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Fig. 2. — Example DEIMOS spectra of Bw-band dropout candidates. The spectra are ordered in increasing /-band magnitudes from 
top to bottom. Solid vertical lines mark the location of Lya while dotted vertical lines mark, from left to right, Si II 1260 A, Ol+Si II 1302 
A, C II 1335 A, Si IV 1394,1403 A, Si II 1527 A, C IV 1548,1551 A in absorption, respectively. 



than that at z = 3.3 (4.1) by ~ 0.08 dex (assuming that 
the UV spectrum has the form f\ oc A -1 " 45 ). 

In order to create average SEDs of galaxies within each 
luminosity bin, we first constructed stacked two dimen- 
sional images for each band. To create clean stacked im- 
ages in each band, unbiased by outliers and sky subtrac- 
tion uncertainties, we adopted the following procedure. 
(1) We first convolved all the images in a given band to 
a common PSF. (2) Then, from these convolved images, 
we extracted l'xl' image cutouts centered on each can- 
didate. (3) Next, in order to obtain an unbiased estimate 
of the sky brightness level, we generated an object mask 
to mask all sources present in each cutout, including the 
central one (IRAF obj masks task was used to detect ob- 
jects above a threshold then several pixels were added to 
grow the mask ensuring all the "light from sources and 
their halos are effectively excluded). We then estimated 
the sky level using the modal value of the unmasked pix- 
els in each frame, and subtracted this value from the 
image (i.e., to create an image with zero sky). (4) Since 
the image cutouts come from survey data with different 
zero points, we then rescaled each cutout to a common 
zero point magnitude. (5) Finally, we constructed a me- 
dian stacked image for each band from all the scaled, 
sky-subtracted cutouts of candidates within a given lu- 
minosity bin. In constructing the median, we used object 
masks to exclude pixels from all objects other than the 
central one. Using the mean instead of the median does 
not change our results significantly. Before the final step 
of stacking, we visually checked all the cutout images. 
Whenever the central source is significantly blended with 
foreground object, we exclude the source from stacking 
in order to avoid contamination to the stacking signal 
by interlopers. The effective number of sources used in 
stacking is reduced by 5 — 10% depending on passbands 
due to confusion. 

This procedure is repeated for each passband 
(UB W RIJHK S and IRAC bands). Because the blend- 
ing and crowding vary greatly with passbands, slight ad- 



justments were made to the IRAF obj masks parameters 
to ensure that all light from the surrounding sources are 
well masked without masking the majority of the pixels 
in each cutout. The latter would result in substantial 
compromise in the effective depth of the stacked image 
as most pixels are masked and thus not used in the stack- 
ing. In Figure 2] we show the stacked images of all six 
subsamples in the passbands where flux was detected. 
The stacked LBT [/-band data yielded only an upper 
limit as no detection was made. The upper limit ranges 
over AB=29.8 - 30.9 for all of our subsamples (see Table 
-)■ 

We measured aperture photometry from each stacked 
image and thus derived the average SED of each sub- 
sample. The photometry was measured in a circu- 
lar aperture 5" diameter,, then aperture-corrected to 
derive total magnitude. For the optical passbands, 
which are convolved to a common Moffat profile PSF 
(MFWHM=1.35", 0=2.5: see IBrown etafl 127501 . and 
for the near-infrared data, also convolved to a common 
PSF (MFWHM=1.35" for H and K s ; MFWHM=1.6" 
for J, 0=2.5; A. Gonzalez et al., in prep.), the aper- 
ture correction factors were estimated analytically. For 
the Spitzer IRAC data, we use the values published in 
lAshbv et al.l ((2009). The aperture correction factors are 
given in Table 1 in units of magnitudes. 

To estimate the impact of the intrinsic spread of the 
range of galaxy luminosities/magnitudes on the mea- 
sured median value for the sample, we carried out Monte 
Carlo simulations. Each time, a random 30% of the sam- 
ple was chosen and the stacking and photometric mea- 
surements were carried out as outlined above. We re- 
peated this procedure 500 times to construct the flux 
distribution of a given subsample. The intrinsic scat- 
ter is estimated from this distribution assuming that the 
random errors follow a normal distribution (i.e., they can 
be scaled from the photometric error of the full stack), 
and is then added in quadrature to the photometric mea- 
surement error to estimate the total scatter in the stacked 
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Fig. 3. — The average spectral energy distributions of the six 
UV subsamples are shown together with the best-fit population 
synthesis model assuming a constant star formation history. The 
SEDs were obtained by convolving the best-fit model with the red- 
shift distribution N(z) to account for redshift uncertainties. The 
B\yRIJHKs and iRAC photometry probe out to the rest-frame 
2 fim at 3.3 < z < 4.3. On right y-axis, we also show the char- 
acteristic UV luminosity L* at 1700 A z ~ 3 ISteidel et alj|1999) . 
corresponding to the /-band magnitudes (at A obs S3 8100 A, the 
third datapoint from left for each SED). The overall shape of the 
SEDs is remarkably similar to one another suggesting the regularity 
of their formation processes (see text). 

photometry (Table 1). 

As shown in Figure [31 the average rest-frame opti- 
cal SEDs of galaxies in different UV luminosity bins are 
strikingly similar to one another in their overall shapes. 
All of them share the flatness in the rest-frame optical 
portion of the spectrum (/„ oc A ), and exhibit a weak 
Balmer break, which is located between the H and K$ 
bands at redshift z = 3.1 — 4.1. Furthermore, it is ev- 
ident that the UV-brighter galaxies are also brighter in 
the rest optical by roughly the same scale factor as in 
the UV. Despite the overall similarities, however, there 
are some measurable differences. In the rest-frame UV 
portion of the SED, the UV slope /3 (= dlogF\/dlogA; 
measured in the observed RIJ bands) is steeper with 
increasing UV luminosity. In the next section, we dis- 
cuss the physical implications of our results using stellar 
population model fitting of the averaged SEDs. 

5. AVERAGE PHYSICAL PROPERTIES 

5.1. Stellar Population Synthesis Analyses 

We determine the average physical properties, specifi- 
cally the stell ar mass, age, and reddening , of our samples 
by fitting the IBruzual fc CharloU ()2003l ) stellar popula- 
tion synthesis models to the stacked photometry. We 
ass ume the Chabr i er (20 03) initial mass function (IMF), 
the iCalzetti et ail (|2000f) form for the effective dust ex- 
tinction law, and consider only solar metallicity popu- 
lations. We model the star formation histories (SFH) 
to be constant, exponentially declining (SFR oc e _t / T ) 
with a wide range of r values (10 Myr - 5 Gyr), or lin- 
early increasing. The stellar age is allowed to vary be- 
tween 50 Myr and the age of the Universe at the red- 
shift z where the fitting is performed (i.e., « 1.66 Gyr at 
z m 3.7). The minimum age of 50 Myr is set based on dy- 
namical timescale arguments when considering the mea- 
surements of half-light r adii and velocity dispersions of 
galaxies at high redshift (|Bouwens et al .1 12004; Erb et all 



2006). While these assumptions for age ranges, IMF, 
metallicity, and extinction may not be accurate descrip- 
tions of the properties of high-redshift galaxies (and may 
even vary with galaxy mass and/or luminosity), our main 
interest lies in comparing the relative average properties 
of galaxies in different UV luminosity bins. Our cur- 
rent study is limited by the nature of our data, which 
is confined to broad band photometry; more details will 
emerge only from deep spectroscopic observations, which 
is now within grasp for the bright candidates within our 
sample. 

We estimate the uncertainties on the derived physical 
properties that result from the finite width of the red- 
shift distribution using a Monte Carlo approach. We 
randomly draw a redshift from the observed N(z) dis- 
tribution (right bottom of Figure [lj, assign it to the ob- 
served average SED, find the best-fit parameters at each 
time, and repeat this procedure 100 times. Although we 
used the observed spectroscopic redshift distribution for 
these analyses, we also simulated the expected redshift 
distribution given our adopted UV color selection crite- 
ria and a large range of theoretical galaxy SEDs. The 
two distributions agree with each other reasonably well 
(Figure [l}, and we have verified that using one or the 
other does not affect our main results. 

5.2. Luminosity, Mass and Age 

Given the narrow redshift distribution, the galaxies 
in our sample range in rest-frame UV luminosity from 
M 1700 ~ -21.43 to -23.30 (see Table 1). Our sam- 
ple therefore consists primarily of L > L* galaxies, and 
greatly increases the number of candidates for the UV 
brightest galaxies. Based on our population synthesis 
analyses, we find that average galaxies in our sample span 
roughly a decade in stellar mass, from log(M*/M0) ps 
10 — 11, and the best-fit population synthesis ages of 
200—400 Myr in the UV bins considered here. 

We summarize our results in Figure [5] when assuming 
the constant SFH model, which facilitates direct com- 
parison with other studies in the literature. However, we 
note that linearly rising SFH models generally provide 
equally good fit to the stacked SEDs. The average stel- 
lar mass of our samples is tightly correlated with the UV 
luminosity, as expected from the homologous SED shapes 
in the different luminosity bins (Figure [3]). Assuming the 
median redshift z — 3.7 to convert the /-band AB mag- 
nitude to the observed UV luminosity, we find that the 
correlation is well described by 

logM*/M = -0.413 M 1700 + 1.367 

(top panel of Figure [5]) . Since the uncertainties in stel- 
lar mass fully account for both random and intrinsic 
scatter in the photometry within the sample, and the 
width of the redshift distribution, the intrinsic correla- 
tion is likely even tighter than measured here. The UV 
luminosity range of the NDWFS sample overlaps with 
the luminous e nd of the sample (about 70 - 80 galax- 
ies) probed by IStark etHI (|2009f ) at M uv < -21.2 at 
similar rcdshifts. Our estimation of the median stellar 
mass at these UV bins is slightly higher than theirs as 
shown in Figure (dashed line; we adopted the stel- 
lar mass when constant star formation is assumed), but 
considering the large scatter and small number of galax- 
ies included in these bins, the two measurements are in 
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TABLE 1 

Average Photometry of Star-Forming Galaxies atz~ 3.7 



Average Photometry 





Am a 


it 1 




=tt.3 

7T° 


it 1 


*5 


Jin 




_ 


19.46 < I b < 23.06 


23.06 < I < 23.66 


23.66 < I < 23.96 


23.96 < J< 24.26 


24.26 < I < 24.46 


24.46 < I < 24.96 


Ml700 c 




-23.30 


-22.51 


-22.18 


-21.89 


-21.63 


-21.43 


N d 




74 


254 


322 


464 


307 


495 


U 




> 29.84 


> 30.51 


> 30.63 


> 30.83 


> 30.61 


> 30.87 


B w 


0.09 


25.97 ±0.18 


26.93 ±0.27 


27.01 ± 0.28 


27.42 ± 0.34 


27.86 ± 0.42 


28.73 ±0.63 


R 


0.09 


23.01 ±0.11 


23.68 ±0.06 


23.94 ± 0.06 


24.17 ±0.07 


24.35 ± 0.08 


24.52 ± 0.08 


I" 


0.09 


22.59 ±0.11 


23.38 ± 0.04 


23.71 ±0.05 


24.00 ± 0.06 


24.26 ±0.07 


24.46 ± 0.07 


J 


0.14 


22.21 ±0.18 


23.21 ±0.08 


23.45 ± 0.08 


23.93 ±0.09 


24.12 ±0.17 


24.38 ±0.19 


H 


0.09 


21.92 ±0.23 


22.96 ± 0.11 


23.31 ± 0.13 


23.63 ±0.15 


23.84 ±0.29 


24.09 ±0.20 


Ks 


0.09 


21.56 ±0.20 


22.75 ±0.22 


22.97 ±0.18 


23.38 ±0.25 


23.47 ±0.52 


23.80 ± 0.40 


3.6/mi 


0.24 


21.70 ±0.24 


22.65 ±0.09 


22.84 ± 0.08 


23.22 ±0.07 


23.47 ±0.11 


23.76 ±0.10 


4.5/im 


0.25 


21.74 ±0.23 


22.60 ±0.11 


22.78 ±0.09 


23.11 ±0.09 


23.48 ±0.14 


23.74 ±0.13 


5.8/im 


0.38 


21.57 ±0.29 


22.53 ±0.32 


22.68 ±0.30 


23.01 ±0.36 


23.44 ± 0.65 


23.61 ±0.65 


8.0/jm 


0.51 


21.67 ±0.30 


22.44 ± 0.26 


22.61 ±0.31 


23.31 ±0.39 


23.78 ± 0.72 


23.71 ±0.66 



a Aperture correction factor applied to the photometry measured in a 5" diameter to obtain the total magnitude listed in the Table. 
b The Kron-like total magnitude in the /-band (SExtractor MAG-AUTO: IKron|[T980l ; IBertin & Arnoutsll 19961) 
c Absolute magnitude at 1700 A assuming redshift z = 3.7 
d The number of the photometric candidates used for stacking 

° /-band photometry measured in a 5" diameter with the aperture correction (see §3) 




Fig. 4. — The stacked images in the U B\y RIJ H K s[3.f>][4.5][5.8][8.0] passbands are shown for all six subsamples (the brightest bin on 
bottom). Each image has a dimension of l'x 1' centered around the photometric sample. The central source is clearly detected in all 
passbands out to IRAC 8.0 £im with the exception of the {/-band. 



qualitative agreement with e ach other. We f urther com- 
pare our results with that of IShaplev et a l. (2005), who 
measured the stellar mass and UV luminosity of mod- 
erately UV-luminous star-forming galaxies at z ~ 2.3. 
While they found no correlation between the UV lumi- 
nosity and stellar mass in contrast to our results, their 
estimation of the median stellar mass is very consistent 
with ours, log M*/M Q = 10.0 ± 0.£3, at M uv 21.0, 

15 This is excluding five galaxies in their sample with very red 
UV to optical color IZ—Kg > 3.5. Similar sources at z ~ 3.7 would 
be bright enough to be detected (> 5<r) in the IRAC [3.6/xm] band. 



which is our faintest bin. Further comparison is not pos- 
sible at this time as only «20 galaxies in their sample 
are bright enough to be in ours. 

Unlike stellar mass, we find no correlation between the 
median age of the stellar population and the UV lu- 
minosity. The best-fit age values (200-400 Myr in all 
cases), suggest that the majority of galaxies are dom- 
inated by young star-forming populations - not just in 

Such sources are unlikely to be common in the galaxies in our 
sample as less than 5-10% of our candidates are formally detected 
in the SDWFS data. 
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their rest-frame UV light, but also in their rest-frame op- 
tical to near-infrared light - that is considerably younger 
than the age of the Univ erse at z ~ 3.7 (~ 1.7 Gyr: 
e.g.. iPapovich et al~ll2001l ). Assuming exponentially de- 
clining SFHs would generally make best-fit ages even 
younger. The young ages are inferred from the weak 
Balmer breaks between the H- and /Cs-band, and are 
thus robust against different SFHs assumed for the SED 
fitting. Finally, our results show that there is a positive 
correlation between reddening and observed UV lumi- 
nosity reflecting the increasingly steeper UV slope with 
luminosity observed in the stacked SEDs (Figure[3]). This 
extends to higher luminosity and redshift the trends ob- 
served in other Lyman break galaxy pop ulations (e.g., 
IPapovich etaUl200H: IBouwens et alj|2009j ). 

5.3. Dust Content 

We confirm the previous observations that more UV- 
luminous star-fo rming galaxies a t high redshift are 
dustier (see, e.g., IReddv et al J 12015 [ShapTev et al.ll2Q0H 
IBouwens et al.1 120101 at z ~ 2, 3, and > 4, respec- 
tively) as shown in Figure [5] (both panels). In addi- 
tion to the reddening parameters derived from the pop- 
ulation synthesis analyses, we can also directly measure 
the UV slope f3 = d\ogF\ / d\og\ from the stacked SEDs 
(RIJ bands: rest 1300 - 26 00 A) after correcting for 
the intergalactic absorption (|Madaulll995[ ) which affects 
the short-wavelength bands. From our UV-brightest 
bin to the least luminous bin, the slope changes from 
/3 = -1.16 ± 0.31 (Sample #1) to -1.37 ± 0.14 (#3) 
to -1.78 ± 0.28 (# 6). Shifting the median redshift 
by ±0.2 does not significantly alter our results (< 10% 
in the estimated UV slope (3). Using the Calzetti ex- 
tinction law an d the correlation between the UV slop e 
and reddening (jMeurer et al.lll999t ICalzetti et al.ll2000D . 
these values correspond to the reddening E(B — V) = 
0.22 ± 0.07, 0.17 ± 0.03, 0.07 ± 0.070. The difference be- 
tween the extinction found by the SED population syn- 
thesis fitting and direct measurements of the UV slope 
is likely due to the fact that the stellar age is allowed 
to vary freely for the former but not for the latter (see 
footnote 15). 

Our measurements are largely complementary to th e 
UV slope determination made by IBouwens et al.l (|2009t) . 
While they also measured the UV slope of UV-luminous 
(L > L*) galaxies, the majority of galaxies in their lumi- 
nous samples comes from the GOODS data (0.1 deg 2 in 
total), and thus the total number of galaxies included in 
these bins is much smaller than ours (by a factor of > 10). 
In Figure [5] (right) we show these results together. The 
two measurements are in quantitative agreement with 
each other while there are slight discrepancies at L w L* . 
This may be in part due to the longer wavelength base- 
line used for our measurements (RIJ; out to 1/im) than 
theirs (i'z) where the UV slope j3 steepens (see, e.g., Fig- 
ure [3] near the rest-frame l//m). Furthermore, inclusion 
of the i?-band in the determination of the UV continuum 

16 A(1600 A) = 2.31(0 - $o) (Eq. 9 in ICalzetti et aLll2uu0t) 
is used for this conversion where they adopted /3q = —2.1 for the 
intrinsic UV sp ectral slope assuming a constant star formation his- 
tory for 1 Gyr (Lcithcrcr et al. 1999). The intrinsic slope does not 
change significantly for the same SFH for 300 Myr, the median age 
found for our sample. 



slope may result in, at least for some galaxies, contam- 
ination of flux from Lya in emission, and as a result, 
cause the UV continuum to appear bluer than it is. In- 
deed, some (> 30%) of our spectroscopic sample have 
been observed with Lya in emission. However, a much 
larger spectroscopic sample is needed to make statistical 
estimates of the fraction of galaxies with Lya in emission 
as well as the distribution of the equivalent widths, and 
thus how this will affect the UV slope of each stack. Fi- 
nally, we note that the selection bias of the photometric 
candidates also contributes to the measurement differ- 
ences; for a fixed /-band (detection band) flux, galaxies 
with bluer R — I colors generally have a higher complete- 
ness, and thus the median UV slope of the sample will 
be bluer (steep e r) tha n the intrinsic one (see Figure 3 of 
IBouwens et al.1 (|2009h . for example). However, correct- 
ing for this effect statistically is not straightforward in 
our case, as we use the stacked photometry for the J- 
band to determine the UV slope. Without knowing the 
intrinsic distribution of I — J colors (which is the quan- 
tity we want to measure), it is difficult to statistically 
correct the UV slope measured from RIJ photometry 
based on the R — I colors alone. However, this explana- 
tion is consistent with the fact that the most discrepant 
measurements are where our selection completeness is 
the lowest (»£*). 

When these effects are properly accounted for, the 
two measurements of the UV continuum slope are in 
good agreement. The overall trend is clear; the average 
UV continuum slope increases with UV luminosity over 
a wide range of luminosities (0.02L* — 6L*). 

6. DISCUSSION 

We have identified the largest photometric sample to 
date of the UV-bright star-forming galaxies at z ~ 
3.7±0.4 by applying a standard color-selection technique 
to the broad-band data from the NDWFS. Our final sam- 
ple consists of 1913 galaxies in the luminosity range of 
L* - 5L*, distributed over an area of 5.3 deg 2 . We have 
divided the sample into 6 bins of UV luminosity (i.e., ob- 
served /-band magnitude) and constructed the average 
spectral energy distributions in each bin. As described 
above, we confirm a strong trend between UV luminosity 
and stellar mass. In this section we discuss the implica- 
tions of these results for the nature of the high-redshift 
star-forming galaxy population. 

6.1. The Mass Dependence of the Star Formation Rate 

The fact that the more UV-luminous galaxies are con- 
siderably more massive than less UV-luminous galaxies 
implies a strong correlation between the instantaneous 
star-formation rates (derived from the rest-frame UV lu- 
minosity) and stellar mass (derived from the population 
synthesis fits, and primarily the mid-infrared photome- 
try). Given the young mean stellar ages (200—400 Myr) 
inferred for the (luminosity weighted) stellar population 
in these galaxies (§4), the tight correlation between the 
current forming rates of stars and existing stellar content 
also implies that the formation of the bulk of their stel- 
lar content is likely associated with the current episode of 
star formation, rather than any previous ones. Further- 
more, we infer that the average star formation histories 
of super-L* galaxies are smoothly varying and long (i.e., 
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at least comparable to the width of the redshift distribu- 
tion N(z), or > 400 Myr) rather than peaky/bursty, as 
the latter would lead to little or no correlation between 
UV luminosity (or SFR) and stellar mass. 

A long, continuous star formation implied by our re- 
sults is apparently at odds with the relatively young 
ages suggested by a weak Balmer break observed in all 
UV bins. Within the conventional "r model", simi- 
larly young ages and the tight SFR-M* relation can be 
achieved only if the formation of galaxies began at the 
same point in time. It is, however, unphysical to imagine 
that the majority of the galaxies (of all luminosities) have 
such a synchronized formation redshift (at Zf ~ 4.2—4.5) 
and star formation histories! To maintain a more phys- 
ical picture, we argue that the average star formation 
history of the UV-bright LBG population needs to be 
one in which the SFR increases with time. Even models 
in which the SFR increase mildly with time would keep 
the age of the stellar population relatively young, since 
the light from the galaxy is always dominated by recently 
formed stars, while maintaining the tight SFR-M* corre- 
lation. We note that similar arguments have been made 
in the literature to explain th e same dilemma found for 
z ~ 2 star- forming galaxies (jMaraston et al.l l2010f ). to 
explain the observed change in t he UV luminosity an d 
stellar mass growth with re dshift dPapovich et alj|2010f> . 
and to model reionization (|Finlator et al.l 20101) . 



6.2. Implications for Star Formation Histories 

The observed relation between the star formation rate 
and stellar mass (^>-M t scaling law, hereafter) holds 
the direct constraint on the rate at which the SFR in- 
creases with time. While a more robust determination 
needs to be based on photometric measurements of in- 
dividual galaxies with spectroscopic redshifts (to accu- 
rately measure both the mean and the scatter of the 
relation), we can make a rough estimate based on our 
stacked photometric results. For simplicity, we assume 
tha t all galaxies in ou r sample lie at z = 3.7 and adopt 
the iKennicutH (|1998l ) calibration^ to convert the me- 
dian /-band magnitude to UV SFRs for each bin, af- 
ter correcting for extinction (derived from the best-fit 
SED model of each UV bin). The logarithmic slope 
of the VP-M* relation ranges over 0.8-1.1 and depends 
sensitively on the adopted extinction parameter. When 
the slope is fixed to unity, the best-fit yields the re- 
lation *(z - 3.7) = (240 ± 30 M yr- 1 )[M»/1O 11 M Q ] 
implying that these galaxies double their stellar mass 
every 290 ± 40 Myo • These results are very simi- 

17 In light of our discussion (§6.1) on rising star formation his- 
tory, it is worth mentioning the uncertainty of this commonly 
used conversion on different SFHs: SFR(M0 year -1 ) = 1.4 X 
10 -28 L„(ergs s _1 Hz -1 ). For a fixed IMF, the conversion fac- 
tor is lower by 10-15% (7%) at population age of 100 (300) Myr for 
star formation histories that rise linearly with time. As for expo- 
nentially declining r models, the factor is higher by less than 5% 
for t > 300 Myr at population ages older than 100 Myr. Consider- 
ing the mild dependence of the conversion factor on population age 
and different SFH models, our estimation of the ty-M* scaling law 
should be insensitive to the range of stellar population parameters 
of galaxies within a given UV bin. 

18 The star-formation rate is the rate of increase in stellar 
mass, or ^ = dM*/dt (modulo mass loss). Hence, a linear re- 
lation Af* = T*!f implies M* (t) oc e*/ T for galaxies to maintain the 
scaling relation with cosmic time. The mass-doubling time is given 
as t In 2. 



lar to that found for star-forming galaxies at z ~ 2 
(|Daddi et all 120071) where the \&-M* relation is well de- 
scribed by #(z ~ 2) = (200 Af o yr^ 1 )(M»/lO 11 M ) - 9 , 
and are somewhat higher than that fou nd for IRAC- 
select ed star-forming galaxies at z ~ 3 fMaedis et alj 
2010) for which the scaling relation of ty(z ~ 3) = 
(330 M yr -1 )(M >t /lO 11 M ) - 91 is found. Our results 
suggest that the specific star formation rate (sSFR), on 
average, remains roughly constant or changes only mildly 
at best with stellar mass (see, also. IDaddi et all 1200 
r Pannella et all 120091: iMaedis et al 



[ Dunne et ail [2 009; 



2010t iKajisawa et al 



20101: IReddv' et al . 2006, 2010). 



While our stacking results provide useful insight into 
how the sSFR changes with stellar mass, especially in 
the luminosity regime where no other measurements are 
currently available, we caution readers that the current 
data in hand is not adequate to constrain the full dis- 
tribution of sSFR at giv en stellar mass bins (see, e.g., 
IReddv et all 120061 120101 based on individual measure- 
ments) as one needs to properly factor in the full red- 
shift distribution and K-correction into the error bud- 
get. Furthermore, to infer the scaling law between sSFR 
and stellar mass, we made an additional assumption that 
all galaxies participate in the tight ^P-M* relation with 
no significant outliers. In reality, it is still possible that 
there may be a small fraction of very massive galaxies 
(say, M* = 10 11 M©) that are not very UV-luminous (say, 
I = 24.5) as the star formation started to fade. While 
such a population would not bias our primary measure- 
ments (i.e., the median stellar mass in a given UV bin) as 
long as their fraction is low, they can lower the median 
specific star formation rate at a fixed stellar mass. 

6.2.1. Simple Prescriptions for the Star Formation History 

In order to develop some intuition regarding the po- 
tential of the \I/-M* relation to constrain the underlying 
star formation history of a galaxy population, we inves- 
tigated a few simple toy models in which the SFR is 
independent of mass, but solely a function of time, and 
which start forming stars at a common redshift. While 
such models naturally result in a constant specific star 
formation rate at any subsequent redshift of observation, 
different models will yield different normalizations at dif- 
ferent redshifts, which can potentially be used to discrim- 
inate between models. In particular, we can compare the 
predictions of the toy models to the observed VP-M* re- 
lationsh ips at z ~ 3.7 (thi s work) and at z ~ 2 (from the 
work of IDaddi et all l200l 120071 IReddv et al.ll2006l ). 

We begin by considering four toy models which share 
the common characteristic that star-formation rates are 
smoothly rising (or constant) with cosmic time: 

1. Exponentially Rising Model (EXP): 
*(t) = *oexp[(t-to)/r] 

2. Power-law Rising Model (POW): 
*(t) - *o(t/t ) a 

3. Constant Star Formation (CSF): 

4. Delayed Star Formation (DLY): 
tf(t) = tf exp [-(t-t )M{t/t ) 
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Fig. 5. — Left: The best-fit parameters for the six UV bins are shown. (Top) The mean stellar mass of galaxies rises steeply with UV 
luminosity (or /-band magnitudes) im plying that galaxie s with higher SFRs are also more massive. We also show the formal best-fit to our 
measurements (solid line) and that of Stark ct al. (2009, dashed line). (Middle) Dust reddening also increases with UV luminosity. The 
best-fit extinctio n values from the SE D fitting (filled circles) and those computed directly from the UV slope (open circles; offset slightly for 
clarity) using the Mcurcr ct al. (1999) correlation are shown. (Bottom) The inferred best-fit ages indicate young ages for all six subsamples 
and show little or no correlation with U V luminosity. R i ght: T he median UV slopes ft at different UV luminosities found from this work 
(filled stars) are compared with those of Bouwcns ct al. (2009) shown in filled circles. For this work, the UV slopes were measured from 
the stacked (median) RIJ photometry when the median redshift 2 = 3.7 was assumed. The RIJ photometry corresponds to the rest-frame 
wavelength range of 1300 — 2600 A at 2 = 3.3 — 4.0. As for the Bouwcns ct al. (2009) points measured from a single UV color, both mean 
(solid error bars) and intrinsic scatter lcr (dashed error bars) are shown. The two measurements are in quantitative agreement with each 
other, and show that the average dust extinction of galaxies increases with UV luminosity over a wide range of Muy = —23 to -17. 



The simplicity of these model parametrizations makes 
them an attractive starting point, b efore consider- 
ing more comple x models (see also Renzini 2009; 
iNoeske et al.l 12007). All models have only one free pa- 
rameter (or, in the case of the CSF model, none), which 
facilitates direct comparison between different models. 
In what follows, we give a brief description of our toy 
models. 

Each galaxy is assumed to have formed at z = Zf and 
followed a given star formation history with cosmic time. 
All models are normalized at t = to to have the star for- 
mation rate \I/(io) = ^0 where time t(z) is defined as 
the time elapsed since the formation redshift Zf. The 
location of the same galaxy on the plane with cos- 

mic time can be followed at each time step or z. Similar 
trajectories can be computed for other galaxies of the 
same SFH but with different normalizations ^o- The \f r - 
M* scaling relation at a fixed redshift is then obtained 
by locating all galaxies of different normalizations at the 
same time step t(z). In Figure [6] we illustrate our cal- 
culation for the power-law (POW) model. The ^-M* 
scaling law obtained with any smoothly rising SFH is 
predicted to have a logarithmic slope of unity, given our 
assumption of synchronized formation redshifts for all 
galaxies. Reproducing a shallower or steeper slope re- 
quires either the formation redshift Zf or SF time scales 
(given as cl or t in our toy models) to have a mass depen - 
dence similar to that discussed by INoeske et al.l (|2007| ). 
For example, requiring galaxies with a higher normaliza- 
tion ^0 to have higher formation redshift would yield a 
shallower slope, which may be in qualitative agreement 
with the hierarchical picture of galaxy formation. We 
defer such modeling to future work as the current mea- 
surements (derived from stacked broad-band data) lack 
the constraining power to such precision. 

For the current modeling, we fix the formation red- 



shift to Zf =8 and the redshifts at which galaxies are 
"observed" to z = 3.7 and z — 2 for comparison with 
the measured relationships. The choice of formation red- 
shift is uncertain, but is unlikely to be critical to the 
current discussion. The age estimates from the popula- 
tion synthesis fitting suggest a minimum formation red- 
shift of z > 4.7, and studies of the global star forma- 
tion rate sh ow evidence of galaxy formation at least to 
z ^ 7 (e.g.. iPapovich et "all 120 101) . Because the elapsed 
time is very small at high redshifts, and since we are 
mainly concerned with the observed star-formation rates 
at z ~ 2 — 3.7, choosing Zf = 10 instead of Zf = 8 makes 
little difference to the results. The inclusion of a range 
of formation redshift and more realistic redshift distri- 
bution N(z) (e.g., such as that shown in Figure [lj will 
indeed introduce scatter to the observed \P-M* scaling 
relation. While we postpone more detailed analyses to 
future work, we note that the amplitude of such scatter is 
likely moderate (< 0.1 dex). Given the narrow observed 
redshift distribution (at least for the z ps 3.7 sample) 
and the fact that we are investigating an ^-M* relations 
derived from stacked data, it is unlikely that these issues 
will significantly affect our results. 

6.2.2. Comparison to Observations 

Following this simple procedure, we investigate how 
each model fares against the observational constraints. 
First, we compute the ^-M* scaling relation at z ~ 3.7 
and determine the best-fit parameter (given as a or r) 
by matching the model predictions to our measurements 
(see Figure [7^ for the time evolution predicted by each 
model). Then using the best-fit parameter, we compute 
the scaling relation at redshift 2 to 6 to follow the redshift 
evolution of each model. As the observational constraints 
(of individual measurements) at other cosmic epochs are 
limited at the present time, we mainly focus on the mea- 
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Fig. 6. — Left: Star formation histories of a power-law model is shown at five different amplitudes as an example when the formation 
redshift Zf = 9 and a = 1.7 are assumed. The stellar mass at redshift z, M*(z), to the zeroth order, is the time integral of the star formation 
rate ^ from the formation redshift Zf to z. The vertical line indicates roughly where the median redshift of the Byy band dropouts lies. 
Right: The — M* scaling relation is shown for the same SFHs. Each of the five parallel lines indicates the time trajectory of a galaxy 
following the SFH shown on left. As the galaxy assembles mass, each trajectory travels from left to right in the figure with cosmic time. 
For the model with the lowest amplitude ty, the redshifts of various cosmic epochs are marked. The 4f — M* scaling relation at a given 
cosmic epoch can be obtained by connecting different trajectories (shown for z = 3.8 as a steeper solid line). 



surements at z ~ 3 . 7 (this work), BzK g alaxies at z ~ 2 
(IDaddi et all [20071: IPannella et all 120091) . and BX /BM 
star- forming galaxies at z ~ 2 (|Reddv et al.l I2006D to 
compare with our models. We note that the scaling 
relatio n from individual measurements by IDaddi et al.l 
(2007) is similar to that meas ured from the stacked ra- 
dio data ()Pannella et alj[2009f > of the same BzK sample. 
In Figure [7]d we show the redshift evolution of all four 
models and best-fit model parameters together with the 
observations at z ~ 3.7 and z ~ 2 (shaded regions). The 
la scatter of 0.25 dex is added for the B\y dropouts to 
account for the redshift distribution and measurement 
errors while for BzKs we use the published value of 0.23 
de£3 (IDaddi et all 120071) . As for BX/BM galaxies, the 
scatter is observed to be ss 0.46 dex, about 70% larger 
than that for BzKs. 

It is evident from Figure 03 that the CSF and DLY 
model predict ~ 0.5 and 0.7 dex decline in the ampli- 
tude of the scaling relation, clearly too steep to 
be supported by the observations at z ~ 2 and 3.7. In 
particular, the CSF model requires that the formation 
redshift is Zf w 4.8 (~500 Myr prior to the observa- 
tions) to reproduce the observed relation at z ~ 3.7F1 
The implication of the relatively low formation redshift 
(zf ps 4.8 for the -B^-band dropouts, and Zf « 2.3 for 
z ~ 2 galaxies) is that the galaxy samples observed at 
different redshift bins cannot be connected to one an- 
other through cosmic time. Such a scenario would not 
only contradict long duty cycles suggested by the tight 

19 The quoted v alue of the semi- interquartile range of 0.16 dex 
ii et al.1120071) corresponds to 0.23 dex in the standard devia- 



tion assuming normal distribution. 

20 Note that the exponentially declining r models will generally 
behave in a similar way to CSF models. However, for a given 
[\E>, M»] value, they require the formation redshift to be even closer 
to the time of observations than CSF models because the SFR was 
higher in the past. 



"■P-M* scaling relation but also would require an unrealis- 
tic picture in which the majority of the galaxies observed 
at z ~ 3.7 would suddenly start forming stars at a very 
short time interval around z = 4.8. 

The EXP model provides a good description of the 
observations at both redshifts when r = 0.50 Gyr is as- 
sumed (consistent with the age constraints). A slight 
change of the formation redshift (zf = 7 or 11 instead of 
8) does not change the timescale more than 50 Myr. The 
POW model also provides a good description at z ~ 3.7 
when a — 1.60 (1.70) and Zf = 9 (11) is assumed. It 
is interesting that the exponent a found by simply re- 
quiring it to reproduce the scali ng relation at z ~ 3.7 
is very similar to that found by iPapovich et all (|2010f ) 
who obtained the exponent a by fiting the redshift evo- 
lution of galaxy samples at a fixed number density at 
z = 3 - 8 (see their Figure 2). Unlike the EXP model, 
the POW model predicts a moderate evolution in the 
amplitude of the scaling relation (« 0.3 dex from z = 3.7 
to 2). While the predicted decline is formally not consis- 
tent with the BzK measurements, the real uncertainties 
(random and systematic) are likely under-r epresented in 
Figure [7t> (shaded area). We note that iReddv et all 
(2006) reported much larger scatter (w 0.4 — 0.5 dex) 
for BX/BM star-forming galaxies at the similar redshift 
range. Furthermore, there are other studies that have 
found similarly large or even larger scatter at z ~ 3 and 
a higher normalization than the Bw dropouts in th e scal- 
ing law (jMannucci et alJl2009b FMagdis et al.ll2010Q . Such 
discrepancies in the literature highlight the existing un- 
certainties and difficulties in measuring such quantities 
at high redshift. 

In summary, our analyses suggest that both expo- 
nentially rising or power-law rising models do better 
than constant SFH model. A more definitive answer 
will have to await better data on samples at multiple 
redshifts. The existing data clearly suffers from different 
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Fig. 7. — Left: The redshift evolution of the ty-M* relation is explored for four S FH models from z = 6 to 2 in comparison with the 
observations at z ~ 3.7 (this work) and z ~ 2 (Daddi ct al. 2007; Rcddy ct al. 2006). The parameters a or r and the formation redshift 
zy are chosen for each model to best reproduce the observed relation at z ~ 3.7. The EXP and POW models with a suitable choice of the 
model parameter may provide the reasonable description of the existing data while the CSF and DLY model do not as they predict a steep 
decline in the scaling relation with cosmic time (see text). Right: Four SFH models that best- fit the observed \f — M* scaling relation at 
z ~ 3.7 are shown. All models are normalized to have \P = 100 MQyr^ 1 at z = 3.7. The time evolution of the SFH (before and after the 
time of observations) and, as a result, the \E> — M, scaling law vary dramatically for different models. Hence, more robust determination 
of the observed ^ — M* scaling relations at different redshifts has the potential to elucidate the star formation histories of high redshift 
galaxy populations and their descendants at later epochs. 

oc exp [t/r] where a or r is positive, unless the galaxies 
are observed in a special time of their formation histo- 
ries. If the onset of star formation in these galaxies is 
not coordinated, the direct implication would be that 
the star formation rate is increasing with time. The rate 
of increase in star-formation rates would be better deter- 
mined by observing similar relations at different cosmic 
epochs as discussed in i j6.2.2l In such a scenario, most 
galaxies are expected to have a fairly long duty cycle until 
the process is terminated or slowed down by a feedback 
mechanism. 

A complementary constraint on the duty cycle of 
star formation can be obtained from clustering mea- 
surements: the measured clustering strength (combined 
with a theoretical model) can provide an estimate of 
the dark halo mass hosting the typical sample galaxy, 
and a comparison between the expected number den- 
sity of these halos and the observed number density 
of the sam ple galaxies yields an estimate of the duty 
cycle (e.g.. iMartini fc Weinberg] I20011 iLee et al.l [2000. 
120091: iMcLure et al.l I2009D . Indeed, the long duty cy- 
cle suggested by our results is supported by the clus- 
tering studies of luminous high-redshift galaxies (e.g., 
lAdelberger et al.l 120051: IQuchi et al] I2005D and the in- 
terpretations in relation to the theoretical predictions 
of dark matter ha lo clustering at the same epoch 
(|Conrov et al.ll2006h . 

In contrast, for much less luminous star- forming galax- 
ies at similar redshifts, the e vidence seems to point to 
the contrary. ILee et afl (|2009f) found that the small-scale 
clustering of the faint UV-selected galaxies at z ~ 4 and 
5 (the deepest data with which the clustering of high- 
redshift galaxies was measured thus far) is much weaker 
than that for dark matter halos, and argued that it may 
be a result of short duty cycles (« 10 — 30% or < 0.4 Gyr 



selection biases, the random (and likely significant) 
photometric errors in the rest-optical photometry of 
high-redshift galaxies, and the uncertainties associated 
with SED modeling assumptions. Nevertheless, it 
is worth noting that comparing the observed ^P-M, 
relations at different redshifts to simple models, such 
as those discussed here, has the potential to elucidate 
the star formation histories of high redshift galaxy 
populations. The key to discriminating between, say, 
SFHs that rise exponentially or as a power-law lies in 
unambiguously determining the overall amplitude of 
the VlZ-M* scaling relation at multiple redshifts. While 
interesting, the prescriptions for star formation histories 
described above are very simplistic and one can easily 
imagine that SFHs of galaxies can be mass dependent. 
However, considering more detailed models is beyond 
the scope of the current work, and unwarranted by the 
data. A future determination of the scatter in the rela- 
tionship (through measurements of individual galaxies, 
as opposed to stacked samples) will bear directly on the 
distribution of star formation history parameters (e.g., 
formation redshifts, star formation rate amplitudes, or 
SFH modes) for the galaxy population. 



6.3. What is the Duty Cycle of Star Formation? 

The UV luminosity of a star-forming galaxy is domi- 
nated by hot stars and therefore measures the instanta- 
neous star formation rate of the system. However, the 
existence of the tight ^-M* relation observed within the 
z w 3.7 ± 0.4 star-forming population suggests that the 
duty cycle of star formation must be significantly longer 
than the time range sampled, i.e., > 400 Myr. Further- 
more, the positive correlation between the two quanti- 
ties - dM*/dt oc M* - yields the solution of M* oc t a or 
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at la level) g overn ing their star formation. Furthermore, 
iStark et all (|2009D . who studied the SED properties of 
similar samples of UV-faint star-forming galaxies, found 
little evolution in ages for galaxies at z ~ 4, 5, and 6. 
They argued that their star formation may not be con- 
tinuous, at least at the constant level, as it would lead 
to older ages, on average, for galaxies at z ~ 4 compared 
to z ~ 6. Alternatively, as they argued, their young ages 
could also be produced by rising SFHs in the same way 
as more UV luminous galaxies. However, s uch a scenario 
would be inconsistent with the ILee et all (|2009l ) cluster- 
ing study as the small-scale clustering of the galaxies 
would closely reflect that of dark halos. 

Taken at face value, this apparent difference be- 
tween the implied duty cycles for star formation in 
UV-luminous and UV-faint galaxies suggests that star- 
formation is a luminosity-dependent process such that 
luminous galaxies rapidly grow in mass by continuously 
forming stars (at an increasing rate) while faint galaxies 
flicker with cosmic time as their UV light turns on 
and off at short intervals (« 200-300 Myr). However, 
it is also true that our understanding of UV-faint 
high-redshift galaxies is more limited than that of more 
luminous counterparts as it is more challenging to obtain 
spectroscopic observations even to identify their red- 
shifts, not to mention more detailed physical properties. 
Further investigations are needed to determine whether 
such differences are due to the observational limitation 
or an intrinsic one. 



6.4. The Descendants of the most UV Luminous 
Galaxies at z « 3.7 

The extinction-corrected star-formation rates of L > 
L* galaxies span 50 - 250 M yr _1 at the epoch of the 
observation (i.e., z w 3.7). The observed relation 
suggests that the star-formation rates will increase with 
time (as oc e*/ T or oc t 1 - 6-1 - 7 ; see Figure [7^). Assuming 
an exponential rise, we estimate that these galaxies will 
double their stellar mass every w 300 Myr until (or un- 
less) the star formation is quenched or slowed by a feed- 
back mechanism. In the absence of quenching, by z = 2 
(i.e., « 1.5 Gyr after z = 3.7), L ~ L* galaxies with 
(M*) w 10 10 M Q will grow into M* « 10 nA M Q while 
2.5L* galaxies will become in excess of 10 12 M Q . Hence, 
the descendants of these UV-luminous star-forming sys- 
tems at z ps 3.7 will likely be galaxies at the top end of 
the mass function by z w 2. 

It is not clear if such growth in the SFR is sustainable 
either by the fuel supply within a galaxy or by the physi- 
cal mechanisms of star formation. More detailed studies 
of individual star-forming galaxies can shed light on the 
gas content and fueling of star formation and the role 
of feedback processes (e.g., AGN). In addition, compar- 
ative statistical studies of the UV luminosity function, 
mass function, and large scale clustering of galaxies at 
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different redshifts (e.g., z s=s 3.7 and 2) can constrain the 
duty cycles of star formation and thereby determine what 
fraction of the UV-brightest galaxies at high redshift will 
have evolved into the most massive galaxies at lower red- 
shift (Lee et al., in prep). We intend to investigate this 
UV-bright sample using both detailed and statistical ap- 
proaches with the goal of elucidating their evolution and 
determining how and when the feedback processes (such 
as AGN feedb ack) will shape the futu re fate of these rare 
galaxies fe.g.. IMarchesini et al.ll2010l) . 
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